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1.  INTRODUCTION  AND  OBJECTIVES 


Sources  cf  coherent  infrared  radiation,  such  as  the  10.6  micro¬ 
meter  laser,  have  become  increasingly  important  in  systems.  Consequently 
there  is  a  need  to  find  improved  techniques  for  beam  control  and  image  or 
signal  processing.  For  example,  there  is  a  liquid  crystal  cell,  but  not 
at  present  an  optical  spatial  lignt  modulator  (SLM)  capable  of  operating 
in  real  time  in  the  3-12  micrometer  wavelength  range. 

Optical  spatial  light  modulations  were  developed  to  transfer 
incoherent  optical  images  onto  a  coherent  light  beam.  Their  use  has  lea 
to  the  meaningful  application  of  real-time  Fouri er-based  optical  signal 
processing  systems.  If  such  a  spatial  light  modulator  were  available  for 
the  infrared  there  would  be  significant  application  there  also.  One 
might  be  artificially  beam  steering  or  focussing  using  the  SLM  as  an 
optically  controlled  diffraction  grating  on  Fresnel  lens.  This  has  been 
demonstrated  by  Christensen  [1]  and  could  have  obvious  applications 
in  an  imaging  array  receiver  used  for  star  system  identification. 

Another  application  might  be  object  identification  by  pattern  recognition 
in  a  10.6  micron  imaging  system.  Such  identi fication  has  been  aptly 
demonstrated  in  the  visible  spectral  range  by  Gara  [2]  using  indoor 
scenes  and  by  Christensen  [3]  using  outdoor  scenes.  Both  workers 
used  a  slow  Hughes  Liquid  Crystal  Light  Valve  type  of  SLM. 

In  yet  another  area,  effort  has  been  put  into  the  development  of 
numerical  optical  computing  using  residue  arithmetic  [4].  The  best 
approach  involves  the  application  of  a  fast  optical  spatial  system  using 
a  cyclic  optical  variable,  the  residue  arithmetic  being  basically  cyclic. 
A  numerical  optical  spatial  computing  system  using  a  bi ref ri ngence-based 
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liquid  crystal  light  valve  has  been  developed  but  is  slow.  A 
seii conductor-based  spatial  light  modulator  using  optical  -araday 
rotation  would  also  have  the  advantages  of  optical  spatial  systems, 
namely  high  information  density  and  parallel  processing,  in  addition  to 
the  cyclic  variable  for  residue  arithmetic  and  the  speed  tnat  other 
systems  lack. 

The  objective  of  this  research  was  to  review  proposed  methods  for 
achieving  optical  spatial  light  nodulation  and  to  make  experimenta1 
studies  of  the  most  promising  approaches. 

2.  REVIEW  OF  SELECTED  OPTICAL  MODULATION  PROCESSES 
2.1  Franz -Kel dysh  modulators 

One  class  of  el ectro-absorpti on  infrared  modulators  is  based  on 
the  Franz -Kel dysh  effect.  This  effect  is  a  variation  of  bandgap  of  a 
semi  conductor  produced  by  means  of  a  high  electric  field.  Devices  in 
practical  use  are  capable  of  an  attenuation  range  of  about  60  dB  and  have 
switching  times  in  the  microsecond  range.  These  devices  are  essentially 
small  optical  waveguides  of  the  order  of  a  millimeter  in  length.  The 
transmission  path  is  a  compound  semiconductor  with  a  bandgap  energy 
slightly  greater  than  the  photon  energy  of  the  radiation  to  be  transmitt¬ 
ed.  When  a  sufficiently  large  electric  field  is  applied  perpendi cul ar  to 
the  guide,  the  bandgap  is  reduced  enough  to  allow  photon  absorption  [7], 
The  physical  embodiment  would  be  an  array  of  small  optical  waveguides 
absorbing  light  at  one  angle  and  radiating  spatial  modulated  light  in 
another  direction.  The  advantages  of  this  system  are  the  existence  of 
the  technology  and  the  wide  dynamic  range.  Di sadvantages  jre  the 
ci rcu i t-impo? jd  1  ini  tat'  ins  on  speed  and  the  geometrical  problems  in 


design  of  arrays.  The  necessity  for  electrical,  rather  than  optical, 
control  of  the  modulation  may  be  a  di sadvantage,  particularly  in  view  of 
the  problem  of  addressing  individual  cells.  For  these  reasons  we  did  net 
investigate  these  devices  further. 

2.2  Sensor-control  1 ed  emission  techniques 

As  the  name  suggests,  such  techniques  refer  to  the  emission  of 
radiation  from  a  source  that  is  controlled  by  light  falling  on  a  sensor. 
In  the  s 1  .op lest  implementation  a  light-emitting  diode  on  one  side  of  an 
opaque  plate  might  be  connected  to  a  phototransistor  on  the  ether  side  of 
the  plate.  Such  an  arrangement  would,  of  course,  require  a  dc  power 
supply.  The  intensity  of  the  emitted  light  could  be  made  higher  than 
that  of  the  controlling  light.  The  emitted  light  could  also  be  of  quite 
different  wavelengtn  from  controlling  light.  If  laser  diodes  were  used 
there  could  be  a  conversion  from  incoherent  to  coherent  light. 

In  order  for  this  sensor-amplifier-emitter  combination  to  be 
useful  in  a  spacial  light  modulator  it  would  have  to  be  produced  in  an 
integrated  array.  Such  an  array  has  been  made  by  Beneking  et  al.  [3] 
using  layers  of  compound  semi  conductors  deposited  by  molecular  beam 
epi taxy . 

Since  we  lack  the  equipment  capability  for  growing  complex  hetero¬ 
structures,  we  made  a  preliminary  study  of  the  feasibility  of 
conventional  integration  methods.  As  a  guide  we  made  several  optical 
logic  gates  as  well  as  amplifiers  using  discrete  components.  The 
approach  is  interesting  for  several  reasons.  The  electronics  of  a  cell 
can  be  complex  while  maintaining  optical  inputs  and  outputs.  Since  the 
output  and  input  are  coupled,  the  addressing  problem  of  the  Franz -Kedysh 
array  is  avoided.  The  dimensions  of  a  cell  might  be  as  small  as  a  few 
microns  and  thus  the  resolution  might  be  adequate  for  many  purposes. 
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Thera  are,  however,  practical  di ff icul ties  in  naming  integrated  arrays  o 
this  type.  They  stem  mainly  from  the  need  for  optical  isolation  between 
elements.  These  difficulties  made  it  clear  that  maxing  sucn  an  a.-ray 
was  beyond  the  scope  of  this  program,  '-leverthe!  ess  toe  fundamental 
advantages  of  electronic  versatility  ir,  an  optical  system  might  justify 
further  work. 

2.3  Free  Carrier  Faraday  rotation  in  semi  conductors. 

Free  carrier  effects  in  semiconductors  can  also  be  used  'or 
optical  modulation.  Tne  basic  idea  is  that  nigh  purity  Gads,  for 
example,  is  quite  transparent  for  infrared  radiation  in  a  band  witn 
wavelength  greater  than  about  1.0  micron.  If  free  carriers  are 
introduced  by  injection  or  radiation,  two  effects  may  be  observed. 
Radiation  at  a  longer  wavelength,  e.g.,  10u ,  will  be  absorbed  in 
proportion  to  tne  free  carrier  concentrati on.  In  addition,  if  a  nagneti 
fie!'  is  applied  in  tne  direction  of  beam  transmission,  there  will  oe  a 
Faraday  rotation  of  the  angle  of  polarization. 

If  an  intense  beam  of  light  of  wavelengtn  0.9-1.0u  is  used  to 
generate  free  carriers,  it  may  be  considered  to  ce  the  control  signa1. 
This  control  signal  can  then  be  used  to  amplitude  modulate  a  1 X  beam 
either  with  or  without  a  magnetic  field.  If  a  magnetic  field  is  present 
it  should  be  possible  to  use  two  coincident  control  beams  to  achieve 
logic  functions. 

A  simplified  theoretical  model  was  presented  in  the  original 
proposal  submitted  on  March  23,  1979.  Expressions  for  the  rotation  per 
unit  length  and  absorption  per  unit  length  were  obtained  from  Moss  [9]. 
Since  then  a  complete  derivation  has  been  obtained  where  classical 
physics  has  been  used  as  a  starting  point,  with  the  results  which  are 
stated  here: 
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where  a  -  absorption  coefficient  per  jr.i*  lengtn 
q  =  electron  charge 


=  wave  lergtn  of  infrarec  rad^..tior 
n  =  Index  of  refraction  of  semi  conductor 
3  =  speed  of  1 1 gnt 
c0  =  permittivity  in  free  air 
un  =  .] act ron  mobi  1  i ty 
yn  =  hole  mob i 1 i ty 
me  =  electron  effective  mass 
nn  =  hole  effective  mass 
5  fr  =  free  carrier  Faraday  rotation  in  RAD/cm 
Bz  =  applied  magnetic  field 
n(x)  =  electron  concentration  in  1/cc 
p(x)  =  hole  concentration  in  1/cc 
In  addition  the  figure  of  merit  can  be  expressed  as  6^/x. 

A  few  comments  can  be  made  about  equations  (1)  tnrough  (3). 

First,  it  should  be  mentioned  that  the  optically  excited  electron  and 
hole  populations  give  rise  to  the  Faraday  effect  are  not  uniform.  Their 
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concentration  profiles  are  usually  exponential  witn  an  associated 
cnaracteri  stic  diffusion  length.  Furtnenore,  .cost  semi  conductors  .nave 
some  type  of  initial  doping.  As  a  result,  n\x)  is  toe  addition  of  toe 
oacxground  level  of  free  electrons  and  toe  optically  excited  ores. 
Secondly,  tne  Faraday  effect  contains  real  and  \megi nary  parts.  The  real 
part  is  tne  expression  for  a  real  rotation  while  tne  i  naginary  part  is 
expression  for  a  real  rotation  wnile  the  imaginary  part  is  an  expression 
for  tne  ellipticity  of  tne  rotated  intrareo  bean,  ""his  elliptic!  t:y  is 
directly  proportional  to  the  absorption  coefficient  which  is  rather  low 
even  for  moderate  levels  of  photo  excitation. 

Several  conditions  are  to  be  net  to  satisfy  tne  above  expressions. 
These  include: 

1)  X  >  ,\g  where  \g  is  the  absorption  edge  wavelength 
of  the  semi  conductor 

2)  <<  (air:n*c)2 /(qB2)^  for  Sz  up  to  10  kGass  th’s 
condition  can  be  easily  net. 


Optical  Generation  oc  Free  Carriers 

The  bul<  rate  of  carrier  generation  at  a  plane  situated  at  a  dist¬ 
ance  x  from  tne  illuminated  surface  of  the  semi  conductor  is  given  by  [10] 

g(x,t)  =  uy\ (1-3.X  )I0(t )s  *ax  carriers/sec-cm3  (4) 

wnere  a  :  absorption  coefficient  in  1/cm 


Y,\  :  quantum  efficiency 

R\  :  reflection  coefficient  of  incident  radiation 
I0(t)  :  number  of  photons/sec^cm2  at  incident  surface 


The  2nd  order  time-dependent  differential  equation  that  describes 
the  motion  of  free  electrons  and  holes  in  a  semiconductor  is  given  by  [11] 
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and  s  =  surface  recombination  velocity  in  cm/sec 
d  =  sampl e  thickness. 

Although  the  above  solution  is  general,  it  does  not  account  for  multiple 
internal  reflections  that  exist  within  the  slab  for  incident  photon 
energies  which  are  just  under  the  semi  conductor' s  band  gap  energy. 

One  can  deduce  from  the  above  formulations  that  due  to  the 
gradients  that  exist  in  the  photo-excitation  of  carriers  along  with 
surface  recombination  phenomena,  the  Faraday  effect  will  be  much  lower 
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tnan  anticipated.  Figures  1  ana  2  snow  the  carrier  distribution  in  a 
250  micron  indium  antimonide  wafer  for  several  input  laser  oower 
densities  (n^  =  1.0)  and  two  different  surface  recombination  velocities. 
Figure  3  shows  the  carrier  distribution  for  a  Ga 1 1 i un-Arseni de  sample. 

Tre  photo-exci ted  carriers  decay  in  GaAs  in  a  shorter  distance  because  of 
their  lower  diffusion  lengths  as  compared  to  InSb. 

'using  equations  (1)  -  (9)  tne  real  part  of  the  Faraday  effect  was 
obtained  for  InSb  as  a  function  of  incident  radiation  power  densities. 
Figure  4  shows  the  angle  of  rotation  for  two  temperatures  in  InSb. 

Figure  5  shows  tne  Faraday  effect  in  GaAs. 

Reduction  in  nobilities  and  lifetimes  due  to  the  high  levels  of 
carrier  concentration  at  high  input  powers  nave  not  oeen  taxen  into 
consideration.  As  a  result,  tne  expected  results  fcr  the  Faraday  effect 
should  be  lower  than  portrayed  by  Figures  4  and  5. 

The  Faraday  effect  is  greater  in  InSb  than  in  GaAs  as  Figures  4 
and  5  snow.  This  is  a  direct  result  from  the  fact  tnat  tne  electron 
effective  mass  is  smaller  in  InSb,  In  addition,  the  carrier  diffusion 
lengths  are  longer  in  InSb;  thus  moderate  levels  of  photo  excited 
carriers  can  be  sustained  for  longer  distances. 

As  we  have  seen,  large  powers  are  required  from  the  control 
radiation  to  obtain  practical  angles  of  rotation.  To  increase  the 
Faraday  effect  the  applied  dc  magnetic  field  may  be  increased  at  the 
expense  of  more  power  applied  to  the  overall  system.  The  proposed 
alternative  is  a  type  of  multi-pass  device  where  the  infrared  radiation 
to  be  rotated  is  forced  to  transverse  the  semiconductor  before  it  exists. 
In  effect,  the  i  parameter  in 
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S  =  V  •  3  •  l 

is  being  increased  where  V,  the  verdet  coefficient,  and  3,  toe  applied  d. 
magnetic  field,  are  uncnanged.  The  multi-pass  device  functions  as  a 
semiconductor  waveguide. 

To  provide  a  practical  picture  of  tne  increased  effect  in  a 
multi-pass  device  where  a  spot  size  of  3mm  x  3mm  in  InSb  at  T  =  iQO'rC  is 
used  as  a  basis  for  comparison. 


(control  beam  focused  in 
an  area  -0.1  cm2  ) 


a  =  /905  microns 
P  =  100  W 
s  =  104  cm/sec 
t  =  . 5  mm 

N0  =  9xl014  atoms /cc 
n ( x )  =  uni  form  1 evel 

The  comparative  Faraday  rotations  for  the  conditions  specified 
above  are  as  fol 1 ows : 


L , I_NGL£  .PA S^S. 


D  of  f 
^  on 


0.13° 

1.4° 


_MULH -PA^SS  __(a. ,H t e rr.a.T.  16°JI 
2.5° 

28.0° 


The  advantage  of  a  multi-pass  configuration  rests  on  tne  total 
increase  in  length  that  the  signal  beam  is  to  travel.  This  length  can  be 
expressed  as. 


l  = 


sin  a 


where  d  is  the  control  beams  spot  size,  and  a  is  the  angle  of  internal 
reflection,  a  must  be  greater  than  the  critical  3ngle  which  is  equal  to 
14.5°  for  InSb. 
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Because  of  the  several  passes  that  the  infrared  beam  takes, 
regions  of  high  carrier  levels,  which  are  usually  close  to  the  surface, 
will  be  tranversed  many  times.  As  a  result,  not  only  is  l  being 
increased,  but  V  which  is  directly  proportional  to  carrier  level  will 
also  be  increased. 

Several  problems  need  still  to  be  studied  for  the  multi-pass 
device:  coupling  devices  and  their  efficiencies,  radiation  angle 
sensitivities,  increased  attenuation  of  the  signal  beam,  and  element  size 
optimization  for  integration  purposes. 

2.4  Optical  Modulation  with  Magnetic  Garnets 

Magnetic  garnets  are  materials  that  can  produce  large  rotations  of 
the  polarization  plane  of  the  light  while  maintaining  an  acceptable 
transparency .  They  appear  very  attractive  for  applications  in  devices 
like  spatial  light  modulators  and  optical  gates  provided  that  the  Faraday 
effect  can  be  modulated  by  external  means,  preferable  optically.  Differ¬ 
ent  attempts  to  achieve  that  are  found  in  the  literature,  like  using 
waveguides  of  garnet  and  an  external  magnetic  field  to  cause  the  switch¬ 
ing  [12]  or  using  a  laser  beam  to  move  and  manipulate  magnetic  bubbles  or 
domains  [13]. 

To  be  useful  as  spatial  light  modulators  or  gates  for  optical 
computing  the  modulation  scheme  must  be  easily  integrated  into  an  array 
and  controlled  in  parallel. 
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Two  of  the  approaches  we  studied  were: 

a)  to  provide  localized  magnetic  fields  to 
control  the  direction  of  magnetization 
of  the  garnet  and  thus  change  the  sign 
of  the  rotation 

b)  to  modulate  some  basic  parameter  of 
the  material  which  affects  the  amplitude 
of  the  rotation. 

a)  The  basic  idea  in  this  study  was  to  produce  a  rectangular 
array  of  magnetic  bubbles  in  a  plate  of  garnet  and  to  remove  or  create 
individual  bubbles  by  optical  means.  The  mechanism  proposed  was  to 
surround  each  bubble  location  with  a  conducting  ring  in  which  a  current 
could  be  used  to  produce  the  necessary  local  magnetic  field.  Two  methods 
were  considered:  One  was  to  place  a  photovoltaic  diode  in  each  ring  and 
the  other  was  to  place  a  photoconductor  in  each  ring.  In  either  case  a 
beam  of  light  on  the  cell  could  be  used  to  switch  the  bubble.  (Fig.  6) 

The  material  most  extensively  used  in  our  study  of  the  effects  of 
local  magnetic  fields  has  been  Gd2.lBi o. 9Pe4. 5 (GaAl ) 0 . 4° 1 2  with  a 
thickness  of  4.7  um  [14],  The  Faraday  rotation  for  this  sample  can  be 
seen  in  Fig.  7  and  its  transmission  in  Fig.  8. 

Since  the  compensation  temperature  is  281°k  the  saturation 
magnetization  constant  at  room  temperature  is  small  and  external  magnetic 
fields  of  17x10”4  Tesla  are  enough  to  reverse  the  Faraday  rotation  in  the 
whole  sample.  A  strong  hysteresis  has  been  observed  in  this  material  and 
once  it  is  oriented  in  one  direction,  it  will  stay  oriented  in  that 
direction  until  a  magnetic  field  in  the  opposite  direction  is  applied. 
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Figure  6  Magnetic  bubble  cell. 
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FIGURE  8  Experimental  Transmission  Coefficient  in  Sample  Garnet. 
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Photol  i  tnographi  cal  ly  defined  metal  -nicroci  rcuits  have  been  used  t 
generate  magnetic  fields  locally.  The  first  idea  was  to  produce  a  local 
magnetic  field  large  enough  to  loca1  ly  reverse  t.ie  direction  of  the  mag¬ 
netization  in  the  garnet  but  the  currents  needed  are  too  large.  What 
appears  to  be  a  much  better  approach  is  to  provide  for  a  nucleation  poin 
in  or  near  the  active  element  defined  in  the  garnet.  This  way  magnetic 
fields  between  10-15  oersted  should  be  enough  to  cause  the  switching  for 
a  given  active  element.  It  nas  been  ooserved  that  evaporating  a  very 
small  area  of  aluminum  near  the  garnet  produces  a  very  effective  nuc '■  sa - 
tion  point.  The  nucleation  points  are  defects  that  ne'p  to  star:.  the 
reversal  of  the  magnetization.  Since  reduction  in  tne  thickness  of  tne 
garnet  reduces  the  energy  per  uni:  area  of  domains  whose  magnet izati or  - 
opposite  to  the  external  magnetic  field  the  etching  of  small  areas  defin 
ed  photol  i  thographical  ly  snould  provide  for  nucleation  po*: ots.  Another 
way  to  produce  a  nucleation  point  i n  garnets  with  composition 
E02.  88i  0.2^0.6^64. 40 1  ?  is  to  evaporate  small  dots  of  gold  and  annea’  tn 
sample  at  500°C  for  15  minutes  in  a  forming  gas  atmosphere,  however, 
tnis  procedure  has  produced  cracxs  in  the  surface  of  the  gannet,  danagi" 
the  s  iple. 

It  thus  appears  possible  to  maxe  a  two  dimensional  array  of 
optical  gates  by  using  a  pattern  of  nucleation  points  on  garnet  in 
register  with  a  pattern  of  semiconductor  current  loops.  An  attempt  was 
made  to  fabricate  the  current  loops  using  a  si  1 i con-on-sapphi re  process. 
Since  a  current  pulse  of  the  order  of  30  mA  is  needed  to  switch  a  garnet 
region  30  urn  in  diameter,  the  design  is  not  simple.  Although  we  were  no 
able  to  demonstrate  a  working  device  with  our  facilities  in  the  time 
available,  the  idea  may  warrent  further  study. 
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was  based  on  the 


Our  investigation  of  tne  second  approach 
well -Known  rotation  of  tne  plane  of  polarization  of  a  linearly  polarized 
bean  of  light  that  occurs  when  it  transverses  a  slab  of  material  and  its 
photon  energy  is  close  to  a  transition  in  which  at  least  one  of  its 
levels,  either  tne  ground  or  the  excited  state,  is  a  multi  pi et  and  it  has 
been  split  (Zeeman  effect)  by  a  magnetic  field  either  internally 
generated  or  externally  applied.  In  addition  to  the  usual  parity,  soin 
and  total  angular  nonent  conservation,  the  conservation  oc  the  total 
magnetic  moment  determines  whether  a  specific  transition  between 
sublevels  is  going  to  be  allowed  or  forbidden  for  each  oc  the  conponents 
of  the  incident  oean,  that  is,  linear  right  circular  and  left  circularly 
polarized.  Tne  simplest  case  arises  when  the  ground  state  is  a  singlet 
(S  state)  and  tne  excited  state  is  a  triplet  (?  state)  (:ig.  3). 

If  we  consider  the  incident  linearly  polarized  beam  as  cr.  no  cf 
right  and  left  circularly  polarized  components,  the  allowed  transitions 
for  both  of  them  correspond  to  different  energy  values  and  thus  their 
index  of  refraction  and  their  absorption  coefficients  will  be  different, 
resulting  in  a  rotation  of  the  plane  of  polarization  for  tne  outgoing 
beam  and  also  some  degree  of  ellipticity.  The  specific  rotation  depends 
among  other  factors  on: 
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1)  Separation  of  the  two  active  sublevels 

2)  Energy  separation  between  the  bean  and  the  center 
of  the  transition 

3)  Linewidth 

4)  Transition  probabi 1 i ties 

5)  Population  of  ground  and  excited  ’eve’s 

Modulation  in  various  degrees  can  be  achieved  by  externally  affecting  ar'y 
of  tne  mentioned  factors  and  each  one  could  be  considered  a  different 
approach.  The  two  approaches  considered  in  this  won\  are: 

a)  Coherent  excitation  which  is  directed  to  fdctor  2 

b)  Incoherent  excitation  which  is  directed  t>  factor  5. 

a)  --  Coherent  excitation  can  be  achieved  when  tne  material  oas  a 
very  narrow  absorption  line  and  it  is  pumped  by  a  circularly  polarized 
aye  laser  beam  whose  monocromaticit,-  is  good  enough  so  that  it  car  be 
tuned  very  close  in  energy  to  the  absorption  line  but  witnout  being 
significantly  absorbed.  Under  these  conditions  tne  electrons  are  exc'ted 
conerently  and  tneir  energy  levels  are  shifted  In  a  process  tnat  has  been 
called  adiabatic  following  [15],  Tnis  produces  large  changes  in  the 
Faraday  rotation  effect  of  a  linearly  polarized  beam  wnose  energy  is  ’n 
the  vicinity  of  tne  absorption  peak.  The  responses  are  very  fast 
(picoseconds)  and  by  adjusting  the  frequencies,  the  controlled  beam  can 
be  made  much  higher  in  power  than  the  controlling  beam,  thus  providing 
gain. 
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d)  --  Incoherent  excitation.  The  previous  scheme  ii;nt  not  oe 
possible  in  the  case  where  the  absorption  line  is  not  very  narrow.  Thus 
the  idea  nere  is  to  pu.mp  tne  system  witn  a  circular1,.,  polarized  bean 
tuned  at  the  center  of  tne  line.  The  absorption  of  the  bean  pumps 
electrons  from  the  ground  level  to  one  of  the  excitea  levels  changing  the 
electron  populations  and  producing  large  changes  in  the  Faraday  rotation 
of  another  bean  tuned  to  a  frequency  close  to  the  absorption  line.  It 
has  been  shown  [15]  that  changes  of  arount  1G£  in  the  population  of  tne 
ground  state  are  expected  to  produce  changes  of  over  2C0t-  in  the  Fa  racay 
rotation  effect  of  some  garnet  materials. 

Tna  material  chosen  for  the  initial  tests  *as  a  magnetic  garnet 
with  composition  Gd?.  lBi  0.  9^4. 5(GaA2  )o.4^12>  tnicx ,  epitaxially 

grown  in  a  non-nagnetic  substrate.  The  reason  r  the  c^o’ce  was  tnat  is 
very  large  Faraday  effect  suggested  a  large  separator,  between  sublevels, 
provided  high  Faraaay  effect  in  regions  of  consid enable  transparency ,  and 
no  external  magnetic  field  was  necessary. 

A  study  of  its  optical  properties  (aosorption,  reflection,  -araday 
effect,  etc.)  showed  tnat  its  rotation  mechanism  was  of  the  di amagneti ■: 
type  (as  in  Fig.  9)  and  the  separation  between  the  absorption  peaxs  of 
both  active  sublevels  was  large.  Nevertneless  the  absorption  peaks  were 
very  wide  due  to  inhomogeneous  broadening  and  it  was  clear  that  coherent 
excitation  would  not  be  significant. 

Incoherent  excitation  was  attempted  in  order  to  partially  change 
the  relative  populations  of  the  ground  and  excited  levels.  This  attempt 
failed  because  the  lifetime  of  the  excited  electrons  was  extremely  short 
( <1  nsc)  possibly  due  to  spin-spin  relaxation  between  the  sublevels  of 
the  excited  state.  This  made  it  impossible  to  populate  one  of  the 
excited  state  sublevels  enough  to  observe  the  desired  effect  despite  tne 
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Figure  10  Optical  light  modulator  based  on 
a  semiconductor  crystal. 


fact  that  the  laser  power  density  was  increased  to  levels  that  damaged 
the  surface  of  the  garnet.  An  attempt  to  find  materials  with  more 
suitable  energy  level  configuration  has  not  yet  been  successful. 


2.5  Modulation  of  optical  absorption  and  refelction  by 
pumping  of  nole  states. 

Infrared  light  todu! at  ion  using  semi  conductors  is  not  new.  Prac¬ 
tical  devices  have  been  made  whe^e  the  excess  carriers  are  generated  by 
minority  carrier  injection  by  a  p-n  junction,  jemaniun  (Ge)  was  used  i.n 
early  modulators  [17],  To  improve  on  its  switching  time  researchers  have 
studied  other  modulators  based  on  direct  bard-gao  ->aterials  sucn  as  Indium 
-Antimonide  [18]  (InSb)  and  Gal  1 ium-Arsenioe  [19]  'GaAs).  Nonetheless, 
these  modulators  required  the  fabrication  of  p-n  junctions  to  generate 
tne  holes  needed  for  light  absorption  near  ten  microns  y;n).  Other 
researchers  have  been  able  to  modulate  light  utilizing  another  light 
source  and  a  host  semiconductor.  Specifically,  10.6pm  radiatian  nodu¬ 
lated  by  a  pulsed  ultraviolet  laser  with  a  wavelength  of  0.3371  Mm 
has  been  demonstrated  in  GaAs  [20].  An  SLM  based  on  this  approach  would 
almost  be  impossible  due  to  its  complex  arrangement.  Germanium  couplers 
are  utilized  to  direct  a  10.6  um  beam  in  and  out  of  a  2.0  millimeter  long 
GaAs  waveguide.  More  recently,  modulation  of  3  pa  light  by  13  :n  lignt  in 
p-type  Germanium  was  shown  [21].  Specifically,  the  transmission  of  the  3  ,, 
laser  beam  is  increased  by  the  saturation  of  states  on  top  of  the  valence 
band  by  the  10  un  light.  A  practical  modulator  may  not  result  using  this 
schene  due  to  the  rather  large  control  beam  energy  densities  {-  0.5 
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joules  per  squared  centimeter).  A  simpler  approach  to  obtain  an  infrared 
SIM  would  be  to  optically  excite  free  electrons  and  holes  using  a  control 
laser  emitting  photons  with  energies  slightly  larger  than  the  semicon¬ 
ductor's  band-gap  energy.  In  turn,  these  free  carriers  absorb  the  10  j  m 
radiation.  It  has  been  shown  that  the  dominant  absorption  mechanism  for 
light  having  a  wavelength  near  10  yin  is  direct  sub-band  transitions  within 
the  valence  band  [22,2 3],  As  a  result,  low  levels  of  optically  excited 
holes  are  sufficient  to  provide  reasonaole  modulation  depths.  Conse¬ 
quently,  required  control  eneryy  densities  are  expected  to  be  low. 

A  simple  representati on  of  an  optical  light  nodul ator  based  on  the 
absorption  of  light  by  optically  excited  free  carriers  is  oepicted  in 
Figure  10.  Speci f i cal ly ,  1.06  ym  light  supplied  by  a  pulsed  laser  is  to 
be  used  to  control  the  intensity  of  10.59  y'i  light  within  a  host  semi¬ 
conductor.  Silicon  is  chosen  for  this  application  because  of  its  low 
absorption  coefficient  at  the  1.06  yn  wavelength.  As  a  result,  a  uniform 
distribution  of  carriers  can  be  generated  by  the  controlling  laser.  This 
in  turn,  maximizes  the  absorption  of  the  10.59  ym  radiation  which  ’s 
directly  porportional  to  the  product  of  the  number  of  free  holes  and  ics 
optical  path  length.  Otner  materials  with  faster  recombination  lifetimes 
(i.e.  direct  band-gap  semiconductors)  or  silicon  with  gold  centers  can 
used  to  increase  the  modulator's  switching  speed. 

In  Section  2.51,  the  theoretical  analysis  describing  the  ’nter- 
action  of  both  electromagnetic  radiations  with  the  host  semi  conductor  is 
presented.  This  theory  is  separated  into  three  parts;  carrier  generation 
by  1.06  ym  light,  carrier  redi stribution  and  decay,  and  finally,  10.59  ym 
light  absorption.  The  resulting  equations  are  derived  from  the  funda¬ 
mental  relations  such  as  Maxwell's  equations  for  the  carrier  generation 
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and  absorption  theory,  and  the  continuity  equations  for  the  carrier 
recombination  theory.  Infrared  absorption  is  describee!  by  tne  exponen¬ 
tial  decay  of  photon  flux  with  distance  into  the  semiconductor.  Finally, 
a  brief  discussion  of  carrier  redi stributi on  parallel  to  the  sample's 
surface  plane  is  presented.  Tnis  treatment  provides  insight  into 
resolution  limitations  in  a  proposed  SLM  based  on  tnis  infrared 
1 ight -control  led  opti cal -modul  a  tor. 

In  Section  3.2  tne  experimental  apparatus  utilized  and  the  met.iod 
of  data  collection  is  described.  The  experimental  configuration  is 
presented  in  which  a  silicon  wafer  is  simultaneously  irradiated  by  the 
controlling  beam  and  the  probing  bean.  The  sample  is  adjusted  for 
Brewster's  angle  of  incidence  of  the  10.59  an  beam  whereas  the  1.06  an 
controlling  beam  enters  the  sample  from  the  opposite  surface  norma’  to 
it.  Also  presented  is  a  comparison  between  theoretical  a^d  experimental 
results.  In  addition,  an  improved  optical  modulator  with  a  faster 
response  time  is  included.  Finally,  the  results  of  preliminary 
resolution  studies  of  a  silicon-based  SL'i  are  shown. 

Finally,  Section  4  contains  a  review  of  important  aspects  of  tne 
work.  In  addition,  recommendations  are  made  for  further  studies  includ¬ 
ing  proper  resolution  studies  for  a  spatial  modulator.  Furthermore,  pro¬ 
posals  for  additional  work  are  made  where  other  semiconducting  materials 
could  replace  silicon  as  the  host  crystal  in  an  optical  spatial  lignt 
modul ator. 

The  mathematical  details  of  the  discussions  are  presented  in  a 


series  of  appendices.  The  redistribution  of  the  photo-excited  carriers  is 
described  by  an  ambipolar  diffusion  equation.  The  derivation  of  the 


time-space  solution  of  such  equation  is  covered  in  Appendix  A.  The 
mobility  and  diffusivity  of  these  carriers  are  dependent  on  several 
parameters  including  carrier  density  level.  Appendix  3  covers  the 
complete  formulation  of  an  equation  describing  nobility  as  function  of 
carrier  density.  Finally,  Appendix  C  groups  all  of  the  computer  programs 
utilized  in  this  work.  Data  acquisition  programs  are  also  included. 


2.5.1  Theoretical  investigation 

Semiconductors  may  be  used  as  host  crystals  where  el ectronagneti c 
radiation  can  interact  with  other  radiation  to  obtain  optical  switching. 
By  stimulating  the  crystal  with  an  incident  light  pulse,  el ectron-nol e 
pairs  are  generated  provided  that  the  photon  energy  is  larger  tnan  the 
semiconductor's  oand-gap  energy.  A  second  electromagnetic  wave  traveling 
within  the  semiconductor  sample  whose  photon  energy  is  much  less  than  the 
band-gap  energy  would  then  be  absorbed.  Several  absorption  mechanisms 
are  present  when  the  probing  beam's  wavelengtn  is  near  ten  microns. 

These  include  phonon  absorption,  free  carrier  absorption,  indirect  free 
hole  transition,  and  direct  intravalence  transition.  As  the  excited  cree 
carriers  return  to  their  equilibrium  state  by  recombining,  the  probing 
electromagetic  radiation  returns  to  its  steady-state  value.  As  a  result, 
optical  switching  is  obtained. 

The  purpose  of  this  section  is  to  develop  theory  that  relates  the 
intensity  of  a  control  beam  pulse  to  the  transmission  and  reflection  of  a 
probing  beam  via  the  electron-hole  pair  recombination  within  a  silicon 
crystal.  In  part  A,  the  generation  of  electron-hole  pairs  as  a  function 
of  the  control  beam's  intensity  and  transmission  is  developed. 


27 


Since  tne  photon  energy  of  a  lld:YAG  laser  is  slightly  greater  than  the 
bana-gap  of  silicon,  it  is  assume J  that  for  every  photon  absorbed,  one 
electron-hole  pair  is  generated.  Furthermore,  since  the  absorption 
coefficient  in  silicon  at  1.06  yn  is  less  than  100  cm'*,  the  control 
bean's  energy  will  not  be  absorbed  within  a  few  microns  of  tne  crystals 
surface. 

In  part  B,  a  solution  to  the  diffusion  equation  is  developed 
for  which  the  electron-hole  distribution  found  in  section  A  is  the 
initial  condition.  A  simple  diffusion  equation  is  obtained  provided  that 
the  excited  electron  arj  hole  concentrations  and  their  respective  fi-st 
and  second  derivatives  are  approximately  equal. 

In  part  C,  an  explanation  is  .mven  describing  tne  absorption  o‘ 

10  urn  light  in  semi  conductors.  A  graphical  representation  of  the  a’ lowed 
optical  transitions  in  p-type  silicon  is  included.  These  absorption 
mechanisms  are  grouped  together  in  an  exponential  form  to  obtain  tne 
transmission  expression  for  tne  infrared  probing  beam. 

Finally  in  part  D,  a  two-dimensional  expression  describing  tne 
diffusion  of  optically  pumped  electron-hole  pairs  along  the  plane 
perpendicular  to  the  direction  of  light  propagation  is  presented.  Since 
two-dimensional  image  processing  can  result  from  such  optically  controll¬ 
ed  infrared  modulators,  phenomena  that  night  limit  spatial  resolution 
is  i nvestigated.  Speci fi cal ly ,  generated  carriers  will  diffuse  a way  from 
tne  spot  of  incidence  with  a  characteristic  lifetime  and  diffusivity 
resulting  in  image  'smearing'.  A  set  of  two-dimensional  graphs  are 
included  for  various  carrier  lifetimes. 
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CARRIER  GENERATION 


r\  « 

El  ectron-noi  e  pairs  are  generated  within  the  silicon  crystal  .vhe1' 
it  is  exposed  to  a  light  pulse  of  wavelength  equal  to  1.06  „n  .  Tni$ 
short  burst  of  energy  is  thus  absorbed  as  it  propagates  through  tne 
crystal.  The  number  of  free  carriers  per  unit  volume  can  be  obtained  as  a 
function  of  the  control  beam's  incident,  reflected,  and  transmitted  irtsnsitie 

Three  homogenous  media  can  be  characterized  by  their  propagation 
factors  <1,  <2,  and  kg  as  shown  in  Figure  11.  Regions  1  and  3  represent 
free  space  while  regi on  3  depicts  the  semi conductor  sample  of  thickness 
d.  Tne  el ectromagneti c  fields  in  media  one  through  tnree  are, 
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The  propagation  factor  for  each  region  is  respect i vely , 


<1  =  A  iu  i 


k2 


;2  (e  2  +  ie  p'^2 


(l.h) 


,i.  n 


k3  =  <1 

The  entire  field  description  is  completed  by  determining  E  ,  E., 

C  c 

and  '<2  from  the  experimental  values  of  the  magnifjdes  of  the  transmission 
coefficient  T0  and  reflection  coefficient  R0,  and  the  incident  energy  of 
the  control  bean.  From  Stratton's  [9]  book  the  trannission  and  reflec¬ 
tion  coefficients  for  an  absorbing  layer  are  respectively, 


T 


o 


!  ^i(<2-<3)d 

(l+zi2)(l+Z23)  i  T  rl 2r Zf  ’1<2lj 


(1.9) 


2  i  k  2d 

r12  +  r23E 

Ro  = 

.  r  r  ^2lk2d 

1  +  r12r23e 

where, 

ujkji 

zji 

y  (,  kj 

1  -  Zjt 

1  +  Zji 

rj* 

(1.11) 

l  =  1,2,3  (1.12) 
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transmission  coefficients  as  a  func  t  ion  of  sample  t  hickness 


and 

H  =  -H  *  in-  u.: 

Tie  computed  values  of  T0  and  R0  cam  be  matcned  to  tneir  experi¬ 
mental  values  by  varying  ;  ^  H  k?.  Tnus  i?>  the  absorption  coeff’cie° 
in  silicon  at  a  wavelength  of  1. 06un  can  be  found  experinenta1 !y.  As  a 
result,  the  experimental  values  of  tne  absolute  values  s  loaned  of  E?  a 
£2  can  be  obtained  with  tne  knowledge  of  the  transmission  coefficient, 
neflection  coefficient,  and  incident  energy  of  the  control  beam.  Ca ra 
must  oe  exercised  when  obtaining  tne  experimental  val  je  for  -t ?  when  jy' 
equations  (1.9)  and  ,1.10).  The  expected  val  fpr  S'licn'-  at  tne 
wavelength  of  interest  (1.06  un)  is  of  tne  order  of  1C  cm‘:.  -s  a 
consequence,  a  small  error  in  the  ;eas  u-c-ment  of  tne  samp’e  tni  cxness 
leads  to  a  large  uncertainty  in  tne  experimental  value  of  .1 "nis  is 
direct  result  of  tne  reflection  and  transmission  expressions  whi  ;n  are 
oscillating  functions  of  tie  tniccness  a.  figure  12  plots  equations 
(1.9)  and  (1.10)  in  which  has  3seri  set  to  10  cr*. 

The  mean  intensity  of  the  light  pulse  is  described  by  tie  '■eal 
part  of  the  complex  Poynting  vector 


Thus,  the  energy  per  unit  time  and  unit  area  coupled  into  the 
semiconductor  at  z  =  0  is, 
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wnich  can  also  be  stated  as. 


Sa  =  Si  (1  -  R0  -  T0)  (1.21) 

where  Rn  and  T0  are  the  reflection  and  transmi  ssi on  coefficients 
respecti vely . 

Equations  (1.16)  and  (1.20)  can  be  equated  to  equations  (1.17)  and 

2 

(1.21)  respectively  to  form  two  simultaneous  equations  where  ; E 2 !  and 
2 

jE^I  can  be  determined.  The  results  are, 
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(1.22) 


l£2i2 


'1-R0)  c2*2*  .  To 


!  2a  ?d  -2ct2d, 
c  -  6  ! 


(1-R0)  £’2a2d  _  tc 


~2  ;  2a  ?d  -2ct  2d, 

2 


(1.23) 


Finally,  the  electron-hole  distribution  per  unit  time  and  unit 
volume  at  t  =  0  needs  to  be  determined.  Taking  the  real  part  of  tne 
divergence  of  the  Poynting  vector  within  the  semiconductor  gives  the 
energy  absorbed  per  unit  time  and  unit  volume  as  a  function  of  z.  Thus, 


;(z)  =  Re  (|  Es  x  Hs*) 


(1.2-i) 


substituting  equations  (1.3)  and  (1.4)  into  (1.24)  and  differentiating, 


we  obtain, 


,  .  a  2®  2  r ,  2  hx2z  ...  2  2a2z 
;  (z)  =  - -  |E?j  e  +  |t?|  e 


(1.25) 


If  we  assume  that  one  electron-hole  pair  is  generated  with  tne 
absorption  of  one  photon  of  wavelength  equal  to  1.06  un  for  low  carrier 
levels,  then  equation  (1.25)  can  be  set  equal  to  -C(z)hcu  where  C ( z )  is 
the  concentrati on  of  the  carriers  per  unit  time  and  unit  volume  and  tv„  is 
the  energy  of  one  photon.  Thus, 


3J^2  ,  , c+  2  -2a2z  A  ( c- ,  2  >2z 


(jl!  2y  2h 


1  !E5i 


I  Eo| 


(1.26) 


where  [ E 1 1 »  ( E 2 1  »  and  <*2  are  determined  from  experimental  data. 
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8. 


CARRIED  REDISTRIBUTION 


Once  the  short  burst  of  energy  passes  through  the  silicon  crystal 
at  time  equal  to  0+,  the  photo-excited  electron-hole  pairs  will  diffuse  to 
less  populated  areas.  In  addition,  the  electrons  and  holes  will  recom¬ 
bine  with  each  other  with  a  characteri stic  lifetime.  This  annihilation 
process  will  reduce  the  total  number  of  pairs  as  time  progresses.  The 
continuity  equations  that  describe  this  process  are, 


a  n 
5  t 


••n 


(n E ) 


on 


'1  0  7  i 

vi.r, ; 


5  p 
;  t 


=  -hp  T  ►‘P 


(?E)  -  Dp 


(1.28) 


for  electrons  and  holes  respectively . 

Equations  (1.27)  ana  (1.28)  can  be  combined  into  one  single  equa¬ 
tion  that  describes  the  notion  of  electrons  and  holes  within  the  semicon¬ 
ductor.  Several  assumptions  have  to  be  made  to  obtain  such  equations. 

The  optically  excited  electron  and  hole  concentrations  n'  and  ?'  are 
assumed  equal  throughout  tne  crystal.  As  a  result,  tne  difference  between 
the  total  electron  and  hole  density,  n  and  o  respectively,  is  very  sma’l. 
In  addition,  their  time  and  space  derivatives  are  assumed  t)  be  equal.  In 
reality,  however,  as  the  electrons  and  holes  redi stribute,  tne  electrons 
tend  to  diffuse  faster  than  the  holes  since  the  electrons  nave  a  larger 
diffusion  coefficient  than  the  holes  (Dn  >  Dp),  As  a  consequence,  an 
electric  field  develops  between  the  electrons  and  holes  as  their  imbalance 
increases.  Thus,  the  hole  diffusion  is  aided  by  this  field  while  tne 
electron  diffusion  is  slowed  down.  The  basic  assumption  states  that 
although  tne  electron  and  hole  diffusion  coefficients  are  different,  the 
electric  field  from  charge  unbal 3nce  prevents  the  excited  carriers' 
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concentrati ons  and  their  respective  derivatives  fron  differing  signifi¬ 
cantly  from  each  other  throughout  the  semi  conductor.  As  a  result,  a  sim¬ 
ple  diffusion  equation  with  its  corresponding  “anbipolar"  mobility  ua> 


lifetime  ta>  and  diffusion  coefficient  Da ,  may  be  used  [25]. 

5y  multiplying  equation  (1.27)  by  upp  and  equation  (1.28)  oy  unn, 
adding  the  resulting  two  equations  and  assuming 


op  =  6  n 

■5t  6 1 

p  =  n 

2p  =  2n 

and  using  Einstein's  .Relationship  D  =  —  to  obtain 

qu 

Dpun  ■  Ony  p 
we  may  finally  write, 

5  n  (  ♦ 

_ Upp  +  unn  =  -  Rnupp-  Rpuon  *  unu-.  vp-n,  E  .  n 

6 1 

-  upDntn^p) 

The  electron  and  hole  recombination  rates  Rn  and  Rp  respecti  vely , 
assumed  to  be  linear.  That  is, 


(1.29; 


(1.31) 


.32) 


(1.33) 

are 


^n  =  , 


(1.34) 


and 


.  £l 


(1.35) 


where  n'  and  p'  are  the  optically  excited  electrons  and  holes.  Thu?  the 
total  carrier  densities  n  and  p  are  equal  to  a  background  level  n0  and  p0 
summed  to  an  excited  level  n'  and  p'  respectively.  That  is,  n  =  n'  *  n0 
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H i crone 

Figure  11  Spatial  distribution  and  time  decay  of  excess  carrie 


A  solution  to  the  above  differential  equation  in  one  dimension  , 
may  be  formed  by  a  linear  superpcsi ti on  of  solutions  of  the  form  [26], 


n'  [z  ,t )  =  l  (\  cos  CKZ  +  SK  sine  C^Z)  e 

k=0 


(1.41) 


where 


=  —  +  DaCv^. 
a 


(1.42) 


The  coefficients  A< ,  B< ,  and  C*  must  be  chosen  to  satisfy  the  initial 
condition  (see  equation  1.26), 


o'  (z,t=0) 


5  2“  2  + .  2  -2a  2z  2  2a  2Z 

— 5 Cp;  £  "  +  :Eo  £  c 

•jo  2h  d  1  d 


and  the  boundary  conditions, 


sn' 


■ '  i 


44) 


where  s  is  the  surface  recombi nati on  velocity  at  the  two  boundaries 
depicted  by  z  =  id.  That  is,  the  thickness  of  the  semiconductor  slab  is 
equal  to  2d. 

A  complete  derivation  of  equations  (1.41)  and  (1.42)  is  presented 
in  Appendix  A.  Figure  13  is  a  plot  of  equation  (1.41)  for  an  arbritary 
input  energy  and  an  anbipolar  lifetime  of  4.0  microseconds.  It  should  ce 
pointed  out  that  the  reduction  of  carrier  density  at  the  sample's  surfaces 
is  due  to  surface  recombination.  Here  a  worst  case  value  of  1000  centi¬ 
meters  per  second  is  chosen  for  a  silicon  wafer  250  micrometers  thir<. 

The  time  increments  depicted  in  Figure  13  are  equal  to  3  microseconds. 

As  a  final  note  in  this  section,  the  nobilities  of  electrons  ano 
holes  are  not  assumed  to  be  constant.  As  a  matter  of  fact  tney  are 
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functions  of  the  number  of  impurity  atoms  within  the  crystal,  lattice 
tanperature,  and  the  excess  carriers  [27].  As  a  result,  the  ambipolar 
diffusion  coefficient  will  vary  as  the  number  of  electron-hole  pairs 
decays.  A  complete  formulation  to  obtain  these  parameters  is  included  in 
Appendix  3. 


c.  infrared  absorption 


In  this  section  a  theory  will  be  presented  which  describes  toe 
absorption  of  light  having  a  wavelength  near  ten  microns.  The  absorptior 
of  radiation  energy  can  be  expressed  as  the  exponential  decay  of  tne 
photon  flux  with  distance  insioe  the  semiconductor.  We  nay  write  [28] 

^out  =  ^i n  £  (1*45) 

where  I \ n  and  Iou^  are  the  input  and  output  intensities  of  tne  ten  micron 
light,  at  tne  absorption  coefficient,  and  D  is  the  sample  thickness.  Note 
that  equation  (1.45)  neglects  the  fact  that  part  of  the  incident  light  is 
reflected  at  the  first  surface  where  it  enters  the  sample,  and  that 
another  part  is  reflected  back  into  the  sample  at  the  second  surface 
./here  the  light  exists.  Equation  (1.45)  can  be  satisfied  by  placing 
the  silicon  wafer  at  the  Brewster's  angle  for  the  10pm  radiation. 

There  are  several  absorption  mechanisms  occuring  in  silicon  cor 
radiation  near  ten  microns.  Light  will  be  absorbed  by  the  creation  of 
phonons  within  the  lattice  of  the  semiconductor  [29].  ’he  effect  of 
lattice  absorption  depends  on  the  wavelengtn  of  the  light  and  tne 
temperature  of  the  material.  In  addition,  absorption  due  to  free 
electrons  is  also  possible  although  the  cross  section  for  this  process  is 
much  smaller  than  for  other  absorption  mechanisms.  Finally,  the  mechanism 
in  which  a  free-hole  absorbs  a  photon  by  making  a  direct  intravalence  band 
transition  is  also  included,  these  direct  transitions  include  heavy-hole 
to  light-hole  band  transition,  heavy-hole  to  solit-of*  band  transition, 
and  light-hole  to  split-off  oand  transition  [30], 


As  a  result,  we  nay  write  the  absorption  coefficient  for  silicon 
for  ten  nitron  radiation  as, 

a  ap  +  0Ih+t  +  ah+s  +  a£-'S  +  3  f  (l.'O) 

where  ap  is  the  absorption  due  to  phonons,  the  following  three  terns  are 
the  absorption  due  to  the  direct  intravalence  band  transitions,  and  of  is 
the  absorption  due  to  free  carriers.  The  three  intravalance  band  tran¬ 
sition  coefficients  can  be  grouped  into  one  absorption  tern  a j.  In 
addition,  if  ap  and  of  represent  the  hole  and  free-carrier  radiation 
absorption  cross  section  respectively ,  we  nay  write, 

=  ctp-ra,jp+ofn  '1  . 4  7) 

where  n  and  p  are  the  electron  and  hole  densities.  As  a  result,  the 
transmission  coefficient  of  the  probing  bean  can  oe  rewritten  as, 

T  •  •<»<*'  *  Jfn)“  (1.48) 

where  static  effects  such  as  lattice  absorption  and  any  hole  contribu¬ 
tion  from  doping  atoms  have  been  subtracted  out. 

The  expression  that  relates  the  free  carrier  absorption  coefficient 
to  the  density  of  electrons  within  a  semi  conductor  results  f^on  classical 
dispersion  theory.  The  absorption  increases  as  the  square  of  the  wave¬ 
length  A  and  is  directly  proportional  to  the  carrier  density  n.  We  may 
write  [28] 

A2Q2 

af  =  - 

4ir 2  CJ  n  pc  0 

where  Q  is  the  electron  charge,  C  is  the  speed  of  light,  nr  is  the 


nn  u  n 


(1.49) 
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refractive  index,  e0  is  the  vacuum  perni  ttivi  ty,  mn  is  the  electron 
mass,  and  yn  1S  the  electron  mobility.  This  absorption  is  due  to  tne 
absorption  of  light  wnen  a  free  electron  is  raised  to  a  new  energy 
state. 


The  direct  intravalence-band  transition  can  be  described  as  t*e 
absorption  of  light  when  a  bound  electron  is  excited  to  a  higher  energy 
state  within  the  three  sub -bands  in  the  valence  band.  Both  energy  ana 
momentum  are  conserved  in  the  intravalance  band  transitions,  4s  a 
result,  only  holes  in  a  narrow  region  of  the  heavy  and  lignt  hole  bands 
will  contribute  in  the  absorption  of  13  um  light  [31],  Tne  absorption 
coefficient  is  governed  by  tne  populations  of  these  sub-bands  which  are 
present  when  an  electron-hole  pair  is  generated  by  1.06  urn  light.  Tne 
hole  states  involved  in  the  optical  transition  are  maintained  close  to 
the  dynamic  value  by  various  scattering  processes  of  which  tne  dominant 
is  hole-phonon  relaxation  [32],  This  scattering  rate  occurs  w’th  a 
sub-picosecond  time  constant. 

The  expression  relating  the  absorption  coefficient  as  a  function 
of  available  holes  on  the  valence  sub-bands  results  from  a  quantum  mech¬ 
anical  treatment  of  the  transitions.  The  resulting  form  for  is 
[31,32], 


4-  Q2  p  t 

nrCm  p  I  +  t  zi>1 


(1* 50) 


where  t  and  p  are  the  relaxation  time  and  density  of  holes  respectively, 
and  to  is  the  angular  frequency  of  tne  incident  radiation. 


In  addition  to  the  aforementi oned  processes,  a  dynamic  lattice 
absorption  is  possible  in  silicon.  Two-phonon  absorption  is  energet- 


dependence  of  modulated  CO,,  laser  beam.  Carrier  lifetime  is  4.0  microseconds 


ically  possible  for  a  wave’ engt n  of  10.6  microns  [32],  Toe  probability 
cf  occurrence  i  s  increased  dramatical  ly  since  phonons  are  generates  wn~i 
a  f ree-el ectron  recombines  with  a  nole  during  the  decay  process, 
furthermore,  indirect  free-hole  transitions  are  possible  where  an  elect¬ 
ron  is  excited  within  the  same  band  with  simultaneous  generation  or 
absorption  of  a  phonon  [34].  This  process  follows  the  same  dependence  as 
the  f  ree-el  ectron  absorption  does  where  a  increases  as  a  function  of  ,\2. 
As  a  consequence  of  these  added  absorption  processes,  toe  hole  cross- 
section  oq  in  silicon  is  expected  to  be  slightly  larcer  tnan  auoted  in 
toe  literature  [31]. 

Equation  (1.43;  is  plotted  as  a  function  of  time  in  -ijjre  14. 

The  electron  and  hole  concentrations  are  determined  using  equation  (l.Ai; 
with  the  same  boundary  conditions  utilized  to  plot  Figure  13.  hone  that 
since  the  level  of  photo-excited  earners  is  non-uniform,  an  iterative 
method  with  a  1  micrometer  space  increment  is  implemented  to  obtain  toe 
lornalized  transmission  of  the  10  u’i  light.  This  formulation  assures 
instantaneous  generation  o'  electron-hole  pairs.  As  a  result,  the  rise¬ 
time  of  tne  attenuated  beam  is  close  to  zero. 


Figures  15a  and  15b  show  the  energy  band  diagram  of  silicon 
when  optically  induced  transitions  are  include  [13,20].  The  1.06un  light 
generates  inter-band  transitions  whereas  the  10.59  uni  light  contributes  to 
intravalence  band  transitions.  Note  that  bound  electrons  can  be  excited 
from  the  top  of  the  three  valence  bands  at  <  =  0.  This  is  due  to  the 
small  energy  separation  between  the  split-off  band  and  the  heavy  or 
light-hole  bands.  As  a  result,  1.164  el ecton-volts  (eV)  are  reunited  to 
excite  an  e  >-ctron  f-om  the  top  of  the  split-off  bana  to  tne  conduction 
band.  A  i.06  urn  photon  has  an  energy  of  1.17 2e V . 
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b.  Intravalence  band  transitions  (10.59  urn) 


Figure  15  Silicon  energy  band  diagram  with  allowed  optical 


transitions 


”gy  band 

[10,20]. 
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According  to  Figure  lb,  direct  transitions  are  possible  only 
between  the  split-off  band  and  either  toe  heavy  or  light  hole  band  for 
10.59  pn  photons.  Tne  thermal  energy  kl  at  roon  temperature  allows  toe 
population  of  holes  on  top  of  the  heavy  and  light  hole  bands  to  redistri¬ 
bute  itself.  As  a  result  transitions  1  and  D  'nave  a  nuc*1  larger  probabil 
i  ty  of  occurence  tnan  A  and  3  for  1C. 59  pn  light. 


Figure  16  Carrier  redistribution  normal  to  the 
z-direction.  Lifetime  is  1.0  micro¬ 
second,  grid  size  is  15  x  15  urn?. 


Figure  17  Carrier  redistribution  normal  to  the 
z-direction.  Lifetime  is  100  nano¬ 
seconds,  grid  size  is  15  x  15  um2. 


D.  CARRIER  REDISTRIBUTION  NORMAL  TO  rHE  DIRECTION  OF  LIGHT  PROPAGATION 


In  addition  to  carrier  decay  in  the  z-airection,  the  optically 
excited  electron-hole  pairs  will  diffuse  outward  from  the  area  of  excita¬ 
tion.  This,  of  course,  leads  to  a  three-dimensi onal  picture  of  carrier 
generation  and  decay  within  the  sample.  Equation  (1.40)  can  be  extended 
into  three  dimensions  resulting  in 
fin' (r)  -n' (r) 


fit 


T  a 


-  Da 


‘d'  vr] 


;i  .so) 


where 


r  =  x  ax  y  ay  +  z  a2 


(1.51) 


A  two-dimensional  solution  for  equation  '1.50)  valid  for  t  >  tp  can 
be  wri tten  as  [11,28]. 


n'(x,y,t)  =  - : -  e  4Dat  :a  (1.52; 

TT  t 

where  N  electron-hole  pairs  are  generated  instantaneously  at  a  point 
x  =  y  =  0  and  time  =  0,  tp  is  the  finite  pu’se  wioth  of  toe  exciting 
laser,  and  n0  represents  the  background  level  of  f~ee  carriers.  The 
solution  is  plotted  in  Figure  16  for  an  ambipolar  lifetime  of  1.0  micro¬ 
second,  laser  pulse  width  of  40  nanoseconds  and  300  nicrojoule  energy. 
According  to  the  figure,  a  time  dependent  'snearing1  of  the  incident  laser 
spot  is  apparent.  As  a  result,  the  resolution  of  an  infrared  spatial 
light  modulator  is  dictated  by  material  lifetime  and  nobility.  A  resolu¬ 
tion  of  4  lines  per  millimeter  is  possible  where  the  grid  size  is  15  x  15 
microns  squared.  Figure  17  demonstrates  considerable  improvement  for  a 
tenfold  decrease  in  carrier  lifetime.  As  a  result,  a  material  with  a  100 
nanosecond  lifetime  might  have  a  worst  case  resolution  in  the  order  of  3 
1 i nes  per  mill imeter. 
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3.  EXPERIMENTAL  ST JD IES 

Following  the  theoretical  analysis  of  "araday  rotation  in  semicon¬ 
ductors  (section  2.3)  it  was  decided  to  use  InSb  for  the  experimental 
studies  rat.ner  than  GaAs.  Samples  of  bulk  material  were  obtained  and 
wafers  were  prepared.  Special  attention  was  given  to  characteri zi ng 
samples.  Data  obtained  included  optical  absorption  coefficient  at  1.06 
mi crometers ,  minority  carrier  lifetime,  surface  recombination  velocity, 
electron  nooility  a.nd  mobility  ratio.  The  measurement  techniques  used 
included  Hall  measurements,  transverse  magneto  resistance  measurements , 
photoconductivity  and  photoel ectromagneti c  measurements. *  Since  the 
program  terminated  before  device  studies  could  be  mane  tie  details  of  fie 
measurements  will  not  be  included  here. 

The  basic  experimental  configuration  where  one  1 aser  controls  the 
intensity  of  another  laser  via  a  semiconductor  crystal  requires  that  the 
two  light  sources  be  coincident  within  the  sample.  Few  ni rrors  and  lenses 
are  used  to  direct  and  focus  both  laser  beams  resulting  in  a  very  simple 
experimental  layout. 

This  chapter  is  organized  in  the  following  way:  in  Part  3.1  a 
description  of  the  lasers  used  and  their  respective  detectors  is  presented 
in  Part  3.2  instrumentation  and  circuitry  utilized  in  the  experiment  are 
described;  in  Part  3.3  the  results  of  this  work  are  presented  and  correlat 
ed  with  theoretical  results.  A  step  by  step  procedure  is  also  included  in 
Part  3.3  for  using  some  experimental  results  to  obtain  sample  constants  us 
i ng  theoretical  formulations;  in  Part  3.4  the  improved  switching  time  of  a 
gold  diffused  sample  is  presented;  in  Part  3.5  a  description  of  the  experi 
mental  layout  to  obtain  a  rough  estimate  of  the  infrared  modulator's 
spatial  resolution  is  shown. 

*  Min,  Wi-sik,  Steady  State  Excess  Carrier  Distribution  in  n-Type  InSo 
under  Illumination  by  Light,  M.S.  Thesis,  Ohio  State  University,  1981. 
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Figure  19  Chopped  C0~  laser  light  waveform. 


3.1  APPARATUS  DESCRIPTION 


Figure  18  shows  the  experimental  layout  where  the  sample,  a  silicon 
wa*er  polished  on  both  sides,  is  shown  at  the  center.  Tne  output  from  the 
pulsed  Nd : YAG  laser,  positioned  at  the  top  of  the  figure  is  expanded  into 
a  two  inch  beam  which  is  directed  onto  the  back  side  of  the  sample  for 
normal  incidence  by  a  mirror.  Neutral  density  filters  are  placed  on  its 
path  to  obtain  several  power  levels.  A  lens  is  utilized  to  insure  tnat  the 
1.06  yn  radiation  completely  overlaps  the  10.59  yri  radiation.  An  iris  in 
place  of  the  sample  is  used  to  align  both  beans.  An  RCA  $-1  vacuum  diode 
positioned  in  front  of  the  sample  monitors  the  control  beam's  incident 
and  transmitted  powers.  Finally,  a  beam  splitter  positioned  at  tne 
output  of  the  Nd:YAG  laser  provides  a  trigger  signal  tnrougn  a  Centroni: 

8P.<  65  detector  with  a  1  nanosecond  risetime. 

An  International  Laser  Systems  (ILS)  ^-switched  Neodymium:  Yttrium 
Aluminum  Garnet  (Nd:YAG)  laser  with  a  wavelength  of  1.06  yn  is  jsed  as 
tne  control  beam.  The  40  nanosecond  pulse  provided  by  the  ILS  laser  has 
an  average  peak  power  of  7.5  < i 1 o  watts.  The  probing  infrared  beam  is 
provided  oy  a  modified  MP3  Technologies  carbon-di oxi de  (CO 2)  laser  opera¬ 
ting  at  a  wavelength  of  10.59  urn.  The  CO2  mixture  (1011  CO,?,  10%  Nitrogen, 
balance  Helium)  is  supplied  to  the  continuous  wave  (Cw )  laser  cavity 
utilizing  an  open  flow  system. 

The  signal  beam,  the  10.59  ym  radiation,  is  provided  by  a  chopped  COg 
laser  positioned  on  the  far  right  in  Figure  19.  Tnis  beam  is  directed  and 
focused  to  a  spot  1  millimeter  in  diameter  on  the  .eft  hand  side  of  the 


55 


sample  with  parabolic  and  flat  surface  mi rrors.  The  parabolic  mi rror  is  of 
relative  large  focal  length  to  avoid  any  significant  focusing  as  the  beam 
propagates  through  the  sample.  The  sample  is  positioned  at  the  Brewster's 
angle  to  remove  any  contributions  to  the  observed  material  effect  from  toe 
reflection  due  to  the  and  second  surfaces.  The  wavelength  of  the  CO 2  laser  is 
continuously  monitored  by  an  infrared  spectrum  analyzer  and  Zi nc-Sel enide 
(ZnSe)  beam  splitter  positioned  at  the  output  of  the  laser.  Finally,  a 
pyroelectric  detector  with  a  69  nanosecond  risetime  Field-Effect  transis¬ 
tor  preamplifier  is  used  to  measure  the  time-dependence  of  toe  modulated 
1 9. 59  um  radi ati on. 

The  Nd:YAG  laser  is  pulsed  at  a  rate  of  1  Hertz  to  obtain  a 
stability  within  10%  from  shot  to  shoe.  A  delay  circuit  is  required  to 
synchronize  the  two  lasers  with  the  proper  time  delay.  As  a  result,  the 
simultaneous  occurrence  of  both  lasers  within  the  sample  is  insured. 

Figure  19  shows  the  chopped  19.59  urn  radiation  which  consists  of  a  train 
of  250  microsecond  pulses  at  a  frequency  of  25  hertz.  This  signal  is 
amplified  10,000  times  through  a  low-pass  filter  and  fed  into  a  delay 
circuit  depicted  in  Figure  20.  Consquently,  the  required  1  Hertz 
trigger  pulse  is  generated  by  the  CO2  laser's  output. 

3.2  DATA  COLLECT  ION  AND  INSTRUMENTATION 

A  Tektronix  digitizer  with  its  controlling  microcomputer  with 
graphics  capabilities  is  used  for  all  data  collection  and  storing.  In 
addition,  waveform  integration  and  differentiation  are  possible  through 
numerical  methods. 

Since  the  pulse  energy  of  the  Nd:YAG  laser  is  highly  reproducible 
from  shot  to  shot,  the  incident  energy  on  the  sample  is  separately 
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TABLE 


INSTRUMENTATION 


LASERS 

Type 

tp(measured) 

Power(measured)  X 

C02  a 

ND: YAG 

250  ys  chopped 

40  ns 

2  W  10.59  -n 

7,5  KVJ  1.06  an 

DETECTORS 

Type 

tr 

Responsivity  9  \ 

pyroel ectric 

S-l  diode 

PIN  diode 

60  ns  FET  pre-amp 

1  50  ps 

1  ns  (5 QT. ) 

0.37  nV/w  0  10.59  yi 

20. 0  mV/w  '3  1.06  y  n 

AMPLIFIERS 

Type 

BW 

Purpose 

Keithley  100X 
Keith! ey  1 OOX 

100  MHz  (5  or.) 

10  KHz 

raise  signal  level 
raise  signal  level  for 
delay  circuit 

DIGITIZES 

Tektronix  7912AD 

Tektronix  4051/4631 

500  MHz  (5(E) 

Micro  computer/copier 
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step  1 

Step  2 

Step  3 

step  4 

step  5 

Step  o 


Incioent.reflection,  and  transmitted  1.06  micron  energies 

Adsorption  coefficients  at  l.u6  microns 

iagnituaes  of  electric  fields  squared  at  1.06  microns 
within  the  semicon&ictor 

Time  and  spatial  distribution  oi  optically  generated 
carriers 

Average  nole  cross  section  at  10. bo  microns 
Modulated  time-aepenaent  CU2  laser  beam 


Figure  21  Sequential  steps  followed  for  data 
collection  and  conputation 
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Fiom'p  2?a  Experimental  curves  for  incident  and  reflected  1.06  urn  light  pulses 
normal  to  a  silicon  wafer. 


function  of  1.06  11m  light  energy. 
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Control  Lasar  Energy  Miorojoulee 

Figure  34  Calculated  absorption  coefficient  in  silicon  as  a  function  of  1.06  um 
light  energy. 


measured  from  the  transmitted  energy  by  simply  removing  the  sample. 

These  energies  are  computed  by  integrating  the  laser  pulse  waveform. 

The  intensity  and  baseline  of  toe  CO 2  laser  beam,  is  continuous’./ 
monitored  by  chopping  the  laser's  output.  Because  of  tne  low  respons- 
ivity  of  the  pyroelectric  detector,  a  Keithley  wideband  anp’ifier  witn  a 
gain  of  100  times  is  utilized.  This  signal  is  fed  into  the  digitizer 
triggered  by  the  PIN  diode  sampling  tne  output  of  the  Nd:YAG  laser.  As 
a  result,  the  digitizer  displays  only  a  10.59  gn  signal  *nen  a  1.06  ,ii 
pul  se  is  present. 

Table  3.1  summarizes  all  pertinent  data  corresponding  to  lasers, 
detectors,  and  instrumentation. 

3.3  PRESENTATION  OF  DATA  AND  CORRELATION  WITH  THEORETICAL  RESULTS 

A  silicon  wafer  polished  on  both  surfaces  wi *h  a  thickness  of  250 
microns  is  used  as  a  sample  in  this  experiment.  The  material  is  p-type 
with  resistivity  of  500  ohm-centimeter  correspond! ng  to  a  doping  density 
of  2x10^  Boron  atoms  per  cubic  centimeter. 

Experimentally,  three  sets  of  data  are  generated.  The  order  of 
events  as  depicted  in  Figure  21  is  as  follows:  Incident,  reflected  and 
transmitted  energies  of  1.06  yn  are  measured  from  which  experimental  re¬ 
flection  and  transmission  coefficient  result.  Figures  22a  and  22b  show 
the  resulting  averaged  Nd:YAG  laser  waveforms  for  one  power  level.  5y 
using  equations  (1.9)  through  (1.13)  a  value  fora,  the  absorption  coef¬ 
ficient  in  silicon  at  1.06  microns  is  obtained.  The  resulting  coef¬ 
ficients  are  shown  in  Figures  23  and  24  as  a  function  of  incilent  ener¬ 
gies.  Note  that  close  to  25%  of  the  incident  radiation  is  reflected.  An 
anti refl ecti ve  coating  provided  by  a  thernally  grown  film  of  silicon 
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three  1.06  mn  light  energies. 


300  micro joule  1  .06  vm  light  pulse.  Lifetime  is  4.0  microsec 


r 


Ground  level 


Figure  29  Experimental  time-expansion  of  modulated  CO;?  laser  light  for  three  levels  of 


1.06  yn  Radiation  Constants 


1  .  774  x  101  5 


8.854  x  10-14 


Fa  rads/c. i 


Ml.u2.u3 


4rr  X  TO'9 


Henries/cn 


1.036  x  10-12 


Farads/c:n 


250.0  x  TO'4 


1.871  x  1  0-19  (1  .1  72e'0 


jou1 es 


2.03  x  105 


3.5549  [27] 


0.2982  (X  s  =  -  ) 


dioxide  should  be  able  to  maximize  the  transmission  of  the  1.06  yn  energy 
at  the  first  boundary  [9,13],  Next,  the  absolute  value  of  the  square  of 
the  electric  fields  of  the  1.06  yn  radiation  within  the  semiconductor  can 
be  found  by  using  equations  (1.22)  and  (1.23).  As  a  result,  the  spatial 
distributions  of  the  optically  generated  carrier  density  can  be  computed 
with  equation  (1.26).  Figure  25  plots  the  carrier  density  as  a  function 
of  distance  and  in-cident  energies  for  time  equal  zero.  Table  3.2  lists 
the  constants  utilized  for  the  1.06  yn  radiation  equations. 

Immediately  after  the  light  pulse  passes  through  the  silicon 
wafer,  the  generated  carriers  will  decay  with  an  effective  lifetime  ra. 
Equation  (1.41)  describes  this  redistribution  of  optically  excited 
carriers  and  is  graphed  in  Figure  26.  The  lifetime  of  the  carriers  can 
easily  be  matched  to  fit  the  temporal  variation  of  the  C02  signal  oeam. 

As  a  result,  material  lifetimes  can  easily  be  measured  in  this  contact¬ 
less  method. 

The  experimental ly -observed  nodulated  time-dependent  CO 2  laser 
beam  is  shown  in  figure  27.  A  modulation  near  50%  is  observable  with  a 
control  laser  energy  of  300  mi crojoul es.  Superimposed  on  the  sane  figure 
is  the  time-derivative  of  the  CO2  laser  beam  waveform.  Clearly,  the 
first  derivative  represents  the  occurence  of  the  Nd:YAG  laser  pulse, 
while  the  second  one  of  opposite  polar'ty  shows  the  lifetime  associated 
with  the  decay  of  carriers.  figure  28  is  a  time  expansion  of  Figure  27. 
It  is  clear  that  the  risetime  of  the  modulated  infrared  light  is 
limited  by  the  pyroelectric  detector's  response  time.  In  addition, 
several  different  slopes  are  observed  in  the  CO2  laser  waveform.  These 
correspond  to  the  mul ti -1 i fetime  decay  present  in  indirect  band-gap 
materials  such  as  silicon.  Figure  29  is  a  time  expansion  of  Figure  27 
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1.877  x  IQ'20  (0.1173) 


4.5  x  10-1?  (equation  1.49) 


jou1 ss 


8.14  x  (experimental) 


4.0  (experimental ) 


j  sec 


12 


wnere  tne  three  waveforms  correspond  to  three  different  incident  energies 
from  the  control  laser.  A  saturation  effect  is  noticed  as  the  pump 
energy  is  increased.  This  behavior  can  be  attributed  to  the  fact  tnat 
only  a  fraction  of  the  optically  generated  holes  contribute  to  the 
absorption  of  the  10.59  pm  radiation.  There  is  narrow  region  in  k 
space  away  from  k  =  0  where  direct  transitions  are  allowed.  If  we  move 
away  from  this  region  and  opposite  to  k  =  0,  the  Maxwellian  distribution 
function  for  the  holes  will  become  small,  minimizing  the  probability  of 
transitions  [33,34].  As  a  result,  as  the  electron-hole  pairs  recombine, 
a  good  number  of  holes  are  still  available  on  top  of  the  valence  band  tor 
13.59  yn  absorption.  These  available  holes  are  continuously  shuffled 
around  the  top  of  the  band  by  thermal  scattering  processes  at  picosecond 
rates.  Table  3.3  lists  constants  pertaining  to  the  calculation  of  tne 
absorption  of  3.59un  light. 

Finally,  Figure  29  also  shows  the  theoretical  results  super¬ 
imposed  over  the  experimental  ones  after  readjusting  for  the  saturation 
effect.  Equations  (i.41)  and  (1.43)  are  used  where  a  lifetime  of  4.0 
microseconds  was  selected  to  fit  all  three  waveforms.  In  addition,  a  hole 
cross  section  of  3.14  x  10*^  centimeters  squared  resulted  from  the  ex¬ 
perimental  curves.  This  value  is  of  the  same  order  of  magnitude  as  tne 
ones  published  in  the  literature  (1. 3-7.0  x  10~-5  cm2)  [31,32]. 

From  obtained  experimental  results,  one  can  conclude  that  modula¬ 
tion  of  10  urn  light  may  be  accomplished  with  low  levels  of  optically 
excited  holes.  This  is  due  to  their  rather  large  cross  section  for 
materials  such  as  silicon  and  germanium.  One  should  point  out  that  only 
a  portion  of  these  holes  contribute  to  the  absorption  of  10  um  light.  As 
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a  consequence,  the  modulation  depth  is  not  only  a  function  of  nole 
density  but  it  also  depends  on  tne  number  of  10  pnotons,  Tnat  is, 
for  larger  signal  beam  power  levels,  larger  control  bean  powers  are 
required  to  obtain  similar  modulation  deptns. 

3.4  SWITCHING  IMPROVEMENTS 

An  optical  switch  based  on  silicon  has  one  serious  drawback:  its 
switching  speed.  Due  to  tne  indirect  nature  of  carrier  decay  via 
recombi nati on  centers,  the  lifetime  of  a  free  electron  or  nole  is  in  tne 
order  of  tens  of  microseconds.  It  is  well  <nown  that  this  lifetime  is 
inversely  proportional  to  the  level  of  recombination  centers  in  an 
indirect  band  gap  material  such  as  silicon  [38,39,40].  Tnis  relation  .nay 
be  written  for 

1 

Yn  =  - —  (5*1) 

°nuthni,l 

where  is  the  concentration  of  tne  centers  per  unit  volume,  is  tne’" 
capture  cross  section  for  electrons,  and  v^r,  is  tne  electron  therma1 
velocity.  Similarly,  the  expression  for  the  hole  lifetime  is, 
inversely  proportional  to  the  level  of  recombination  centers  in  an 

1 

t,  = - —  (3.2) 

°pvthp  ^2 

Note  that  Nj  and  Nj  are  the  densities  of  the  two  trapping  levels  assoc¬ 
iated  with  each  gold  atom  [39],  Nj.  is  associated  with  an  acceptor  level 
located  about  0.35  eV  above  the  valence  band  and  N?  is  the  donor  level 
located  about  0.54  eV  below  the  conduction  band.  Tnis  clearly  suggests 
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10.59  um  light 


Figure  30  Experimental  configuration  for  a  gold- 
doped  sample  with  a  metallized  back. 
Sample  is  at  the  Brewster  angle  with 
respect  to  the  10.59  um  light. 
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mat  to  improve  on  the  time  response  of  an  optical  modulator  based  on  p- 
type  silicon,  gold  atoms  should  be  introduced  into  the  crystal  lattice. 

Gold  is  introduced  into  silicon  in  a  controlled  fashion  by 
solid-state  diffusion.  A  chin  film  of  high  purity  gold  (=2Q0A)  is 
evaporated  on  a  silicon  sample  from  the  sane  batch  of  .vafers  used  for 
part  C.  The  silicon  wafer  is  inserted  in  a  Lindberg  diffusion  furnace 
set  at  900C .  (- 3  x  1  01  5  ATQMS/cc)  [41].  In  addi  ti on,  a  thin  film  of 
aluminum  and  gold  is  evaporated  on  one  side  of  the  polished  silicon 
wafer. 

The  metal  backing  on  the  sample  provides  a  two -fold  improvement 
over  the  previous  experimental  set-up.  The  optical  path  of  both  probing 
and  controlling  beams  are  increased.  As  a  result,  an  improvement  in  the 
attenuation  of  the  CO2  laser  probing  beam  is  expected.  Figure  30 
illustrates  the  positioning  of  both  laser  beams  with  respect  to  the 
sample  in  which  the  signal's  incidence  orientation  is  maintained  at  the 
Brewster  angle.  Figure  31  shows  the  time-dependent  waveform  subjected 
to  300  microjoules  of  incident  energy  from  the  Nd:YAG  control  laser. 

Also  superimposed  on  Figure  31  is  the  resulting  theoretical  waveform 
where  a  lifetime  of  0.55  microseconds  provides  the  best  fit.  Thus,  an 
improvement  of  near  one  order  of  magnitude  over  the  previous  lifetime 
(4.0  microseconds)  is  obtained  by  introducing  gold  atoms  into  the  silicon 
1 attice. 

The  number  of  gold  recombination  centers  can  be  increased  by 


lengthening  the  diffusion  time  and  increasing  its  temperature.  As  a 
consequence,  much  faster  switching  times  may  result  (near  1  nanosecond). 
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Figure  n  Spatial  distribution  of  both  10.59  uni  and  1.06  pm  radiation  at  the  sample  surface. 


3.5 


RESOLUTION  ESTIMATE 


A  complex  experimental  set-up  consisting  of  array  detectors  for  Potn 
wavelengths  and  imaging  of  a  standard  resolution  pattern  on  the  sample's 
surface  is  bypassed  in  favor  of  a  simpler  experimental  lay-out  as  shown  in 
Figure  32.  The  basic  premise  is  that  this  will  provide  preliminary 
resolution  data  without  the  need  for  an  involved  experiment. 

The  basic  idea  as  portrayed  in  Figure  32  is  to  focus  the  control 
laser's  beam  with  a  cylindrical  lens  on  the  sample's  surface  depicting  a 
narrow  line.  Then  by  scanning  a  well  focused  probing  beam  in  and  out  of 
the  carrier  generation  region  a  plot  of  the  percent  modulation  as  a 
function  of  space  is  produced.  Clearly,  the  CO2  laser's  spot  size  must  be 
equal  or  less  than  the  Nd:YAG  laser's  focused  line.  As  a  result,  the 
expected  spatial  modulation  follows  the  control  laser's  spatial  intensity. 

The  spatial  distributions  for  both  probing  and  controlling  beams  are 
shown  in  Figure  33.  These  are  obtained  by  scanning  a  pyroelectric  and  a 
PIN  diode  across  expanded  10.59  pm  and  1.06  pn  beams  respectively.  The 
actual  spatial  dimensions  result  from  burnt  patterns  made  on  polaroid 
film.  It  should  be  pointed  out  that  this  method  is  not  reliable  in  measur¬ 
ing  laser  spot  sizes  and  it  is  only  used  to  insure  that  the  10.59  pm  spot 
size  is  less  than  the  1.06  pr.i  light  at  the  sample's  position.  The  norma¬ 
lized  modulation  of  the  CO2  laser  probing  beam  and  the  spatial  distribution 
of  the  Nd:YAG  control  beam  is  plotted  in  Figure  34  as  a  function  of  dis¬ 
tance.  Clearly,  the  modulated  10.59  pm  light  follows  the  shape  of  tne  1.06 
pm  light  indicating  a  direct  transfer  of  spatial  information  from  one  light 
source  to  another.  Unfortunately ,  the  spatial  resolution  is  only  about  2 
lines  per  millimeter.  This,  of  course,  is  limited  by  the  control  laser's 
optics  and  not  the  sample's  characteristics. 
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II 


4. 


CONCLUSIONS 


Several  conclusions  may  be  drawn  from  the  hole  punping  experi¬ 
ments.  To  begin  with,  the  feasibility  of  modulation  of  one  light  source 
by  another  via  the  optical  generation  of  carriers  in  a  semiconductor  has 
been  demonstrated.  This  result  provides  the  basis  for  a  spatial  light 
modulator  for  the  10  u pi  wavelength  region.  Close  to  50%  modulation  was 
achieved  with  relatively  low  energy  density  (6.73  mJcm-’)  originating 
from  the  control  laser.  This  suggests  that  other  light  sources  of 
reasonable  power  output  could  be  used  as  control  sources  provided  tneir 
photon  energies  were  slightly  larger  than  the  semi  conductors '  band-gap 
energy . 

A  saturation  effect  was  noticed  in  the  time-dependent  signal  wave¬ 
form.  This  phenomenon  is  due  to  the  selective  participation  of  free 
holes  in  the  absorption  of  10.59  yn  radiation.  That  is,  only  the  portion 
of  the  free  holes  near  the  top  of  the  valence  band  which  are  located  in 
the  vicinity  of  the  10.59  pn  light  energy  contribute  to  direct  intra¬ 
valence  band  transitions.  As  a  result,  the  attenuation  of  the  CO2  laser 
beam  was  sustained  until  the  non-participating  holes  on  top  the  valence 
band  were  filled  with  the  decaying  free  electrons. 

A  theoretical  fit  to  experimental  data  provided  an  effective  life¬ 
time  of  4.0  microseconds  for  the  high-resistivity  material  used.  By  the 
introduction  of  additional  recombination  centers  into  the  lattice,  the 
carrier  lifetime  was  reduced  to  0.55  microseconds.  This  in  turn, 
improved  the  switching  time  of  the  light  modulator  from  15  microseconds 
down  to  near  1  microsecond  for  a  modulation  depth  of  50%.  Additional 
gold  centers  should  further  decrease  tne  switching  time  into  the  nano¬ 
second  region.  As  a  result,  wideband  optical  modulation  may  be  attained 
with  gold-doped  silicon. 
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According  to  the  simple  resolution  experiments,  a  resolution  of  2 
lines  per  millimeter  can  be  obtained  using  the  4.0  microsecond  material. 
Because  of  the  inability  to  provide  a  well  focused  control  beam  to  the 
sample,  this  resolution  figure  does  not  represent  the  sample's  maximun 
resolving  power.  Further  experiments  are  necessary. 

An  inherent  feature  of  a  silicon-based  optical  light  modulator  is 
its  simple  structure.  An  anti refl ecti on  coating  for  the  control  beam 
could  improve  device  performance  by  reducing  control  beam  energy  require¬ 
ments.  A  thermally  grown  si  1  icon-dioxide  film  could  be  used  for  the 
Nd:YAG  laser  [28],  Simple  fabrication  requirements  are  an  attractive 
featjre  for  spatial  light  modulator  implementation. 

Several  topics  that  merit  further  investigation  have  become  appar¬ 
ent  during  the  experimental  phase  of  this  work.  Lifetime  studies  in  sil¬ 
icon  are  possible  with  this  contactless  method.  Previous  experimental 
configurations  have  relied  on  fabrication  of  p-n  junctions  ta  observe  the 
decay  of  the  generated  carriers.  More  recently,  contactless  methods 
have  been  used  in  which  a  10  Gigahertz  nicrowave  is  utilized  as  a 
probing  signal  [35,36].  Because  all  the  semiconductor  samples  need  to  be 
mounted  within  a  waveguide,  these  schemes  fend  to  be  cumbersome.  It  is 
believed  that  the  configuration  used  in  this  work  provides  ease  of  imp1  la¬ 
mentation  for  lifetime  studies. 

In  addition  to  contact-less  lifetime  studies,  work  should  be 
directed  into  identifying  and  separately  measuring  the  several  absorption 
mechanisms  occuring  in  silicon  near  10  microns.  Futhermore,  competing 
mechanisms  associated  with  the  optical  generation  of  el ectron-hol e  pairs 
and  their  decay  via  recombination  centers  could  be  studied.  Finally, 
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Figure  35  Optical  switch  configuration  based  on  a 
semiconductor  crystal.  10  urn  light  con¬ 
trolled  by  1.06  ym  light. 


Figure  36  Optical  spatial  light  modulator  operating  near  10 

light.  Germanium  film  deposited  on  a  GaAs  substrate 


valence  band  structure  could  be  verified  with  a  tunable  CO2  laser  as  the 
probing  signal.  The  basic  premise  in  tnese  suggestions  is  that  semi¬ 
conductor  processes  can  be  studied  in  a  non-destructive  manner  with 
electromagnetic  radiation. 

A  more  detailea  resolution  study  on  such  an  infrared  SLM  is 
required  to  provide  additional  feasability  information.  Specifical  ly ,  a 
standard  resolution  pattern  should  be  imaged  on  the  surface  of  the  crys¬ 
tal  with  the  help  of  a  lens.  In  addition,  the  control  laser's  beam 
should  be  collimated  and  spatially  filtered  to  obtain  a  clean  wavefront. 
Furthermore,  the  CO2  laser's  probing  beam  should  also  be  collimated  and 
filtered.  Finally,  array  detectors  must  be  employed  to  cet  an  accurate 
spatial  signal  from  such  an  infrared  SLM. 

Several  applications  are  possible  with  a  semiconductor-based 
optical  switch.  The  simplest  one  would  consist  of  a  one-dimensional 
optical  modulator  in  which  a  high-powered  CW  CO 2  laser  is  controlled  by 
a  pulsed  Nd:YAG  laser.  Large  bandwidths  are  possible  if  proper 
lifetime-reduction  techniques  are  employed  in  silicon.  Figure  35  snows 
a  probable  arrangement  where  multiple  passes  of  the  controlling  beam  may 
reduce  its  energy  density  requi rement. 

A  more  complex  impl ementat ion  of  a  spatial  light  modulator  is 
depicted  in  Figure  36.  In  this  particular  instance,  a  10  to  20  pm  film 
of  germanium  (Ge)  grown  or  deposited  on  a  galliun  arsenide  (GaAs)  sub¬ 
strate  could  provide  larger  modulation  depths.  This  results  f ron  the 
fairly  large  absorption  coefficient  of  Ge  at  1.06  yn.  Optical  switching 
is  still  based  in  the  optical  generation  of  free  carriers  which  would 
absorb  radiation  near  10  pm  by  direct  intravalence  band  transitions.  The 
GaAs  crystal  is  utilized  only  as  a  substrate  due  to  its  transparency  at 
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Doth  wavelengths  of  interest.  In  addition,  Ge  can  easily  be  grown  o.er 
such  a  substrate. 
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APPEND  I  (  A 


CARRIER  R£D ISTRISUT IJTJ  [ane  dimension) 

Tne  second  order  differencial  equation  that  describes  toe  deca., 
electron-hole  pairs  nay  be  stated  in  its  ambi, solar  form  as. 


2r 


6  t 


where  n‘  represents  the  number  of  electron -hole  pairs  pen  cubic  centi¬ 
meter  and  :a  and  Da  represent  their  arbi polar  lifetime  in  seconds  ana 
tneir  Jiffusivity  in  centimeter  squared  per  second  respectively. 

If  tne  number  of  carriers  is  expressed  as  tne  product  of  two  in 
pendent  fjnctions, 

n'(z,t)  =  2  (z)  T  (t )  ^ 

we  obtain, 

I  il  =  _  ,  n  1  n  -Z 

T  £  t  :a  "a  Z  3  z  ^ 


by  direct  substitution  of  equation  (A. 2)  into  (A.1).  As  a  result,  we 
write, 


and 


1  62 


-C2 


1  6  T 

T  5 1  "  V 
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wnose  individual  solutions  lead  to  the  overall  expression  for  n'(z,t). 


n'(z,t)  *  'A  cos(cz)  +  B  si.n(cz))e 


•v  l 


(A. 5; 


where 


+  D,C2 


(A.7) 


A  linear  superposition  of  solutions  of  the  form  of  equation  'A. 5) 


sjch  as, 


n'(z.t)  =  I  'A*  cos(CKz )  +  3k  sin(Ckz  ),e  "V!< '  (A.  3) 

k=0 


and 


vk  3  ~  +  °ack2 
T  a 


(A. 9; 


will  satisfy  the  differential  equation  where  the  values  of  AK ,  3„;,  and  C« 
may  be  chosen  so  that  the  boundary  conditions  are  satisfied  by  the 
super-position. 

The  surface  boundary  conditions  are, 


=  0 


a 


z  =  2d 


sn"  (=d,t) 


{A. IQ) 


where  s  is  the  surface  recombi  nati on  velocity  in  centimeters  per  second 
and  2d  is  the  thickness  of  the  semiconductor  slab.  Substituting  a  solu¬ 
tion  of  the  form  of  equation  (A. 6)  into  (A. 10)  we  find  that  in  order  to 
satisfy  equation  (A. 10)  we  need  to  solve  the  transcendental  equation, 

Da 

ctn(Cxd)  =  —  (Ckd).  (A. 11) 
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'»ote  that  for  a  very  good  surface  (s  =0),  C^d  =  where  <  -  0,1  ,2,3, 
and  for  a  lapped  surface  (s  =  >» )  C*d  =  (2k  +  1)  u  / 2  where  k  = 
n  ?  o  ^ 

y  Lm  f  'J  m  m  • 


The  thi-d  boundary  condition  provides  us  with  the  carrier  distri¬ 
bution  as  a  function  of  z  for  time  equal  zero,  (refer  to  equation 
(1.26)). 


n' (z,0) 


-2a2(z-nj)  _  2a2(z+d) 


n  e 


(A. 12) 


.»nere, 


22^2 
H  r  ~~o — 

■x  Hi  n 

By  setting  equation  (A. 9)  equal  to  (A. 12)  we  may  determine  the 
coefficients  A*  and  8i<.  >le  write, 


(A. 13) 


I  i Ax  cos(Ckz)  +  3k  sin(Ckz)j  =  n'(z.O) 

K=0 


(A. 14) 


By  multiplying  both  sides  of  equation  (A. 14)  by  cos(C£z),  integrating 
over  the  entire  thickness  of  the  slab,  and  using  the  orthogonality 
principal , 

d/2 

/  cos(C<z)  cos(C£z)dz  =  0  (A. 15) 

-d/2 

for  k  *  l,  we  obtain  an  expression  for  A^, 
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2Ck 

Ak  ,  _ - 

4 <*2  +  2Ck6  t  sin(2C|<d) 

x  f  C|<sin(Clcd)  {n+  (l-*e  401 2<1)  +n 

+  2a2  cos(Ckd)  {n+  (l-e"4a2d)  0-^)1]  •  (A. 16) 

Similarly,  to  obtain  Bi<,  we  Multiply  both  sides  of  (A. 14)  by  sin  (k $z), 
integrate,  and  use  the  orthogonality  principal.  As  a  result,  we  ha/e, 

1  1 

Bk  =  d  4a |  Ck2 

x  [  C«;  cos(Ct<d)  {n+  (1 -s  ^ 2d)  +  n '  (l-e401^)} 

-2a2Sin(Ci<d)  {n+(l-tc^a2d)-n*(l-e<,a2d)}j  (A. 17) 
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APPENDIX  S 


MOBILITY  DEPENDENCE  ON  CARRIES  LEVEL 

The  /nobility  of  an  electron  or  hole  within  a  semiconductor  is  a 
function  of  three  components:  the  lattice  nobility,  the  impurity  scat¬ 
tering  nobility,  and  the  hole-electron  scattering  mobility  [26].  Tne 
coulomb  interaction  of  the  impurity  scattering  component  and  the  hole- 
electron  scattering  component  are  combined  in  a  reciprocal  law  to  give  a 
mobi 1 i ty  uc, 


1  1  1 

—  =  —  +  -  (3.1) 

he  hi  h  ne 

where  p  i  and  uhe  are  impurity  arid  hole-electron  scattering 
mobi 1 i ties. 

The  effective  mobility  is  obtained  by  combining  uc  and  the  lattice 
mobility  ul  using  a  reciprocal  law, 

1  111 

h  F  h  c  h  L 

where  the  statistical  factor  F  is  defined  as, 

0.  0954  +  0.473Z  -  0.333z2 
F  =  0.0956  +  0.8662  -  0.776 z2 


(3.2) 


(3.3) 
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and. 


^  L 

z  = - . 

UL  + 


(3.4) 


The  expressions  for  pj,  p^e>  and  ul  i n  cm2/volt-sec  are  as 
fol lows: 


w  I 


2.45  x  1321 


tL  1.5 
300 


N  .  In 


1.41  x  1020 

1  +  - 

n+p 


JL  2 

300 


(B.5) 


p  hen 


2.45  x  1021 


300 


H  p  in 


1.41  X1020 

+  - 

n+p 


Ik  2 

300 


(3.6) 


P 

yhep  =  Mhen  “ 

tl  -2*4 

ULP  =  480 


T  -2.5 

ULP  -  480  - 


(3.7) 

(3.8) 

(3.9) 


where  N  is  the  doping  density,  n  and  p  are  the  free  electron  and  hole 
density,  T|_  is  the  lattice  temperature  in  degrees  Kelvin,  and  the  small  n 
and  p  subscripts  refer  to  electron  and  hole  mobilities  respectively. 
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Electron 


According  to  equation  (1.39)  the  anbipolar  diffusivity  may  be 


expressed  as  a  function  of  the  free  carrier  densities  and  their  respec¬ 
tive  di f fusi vi ties. 


p  +  n 

Da  =  -  .  (3.10) 

£_  ■  £ 

Dp  "  Dn 


Figure  37  shows  the  ambipolar  diffusivity  as  a  function  of  free-carrier 
density  level. 
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APPENDIX  C 


COMPUTER  PROGRAM  DESCRIPTION 

Several  computer  programs  are  utilized  in  this  work  to  obtain 
several  coefficients  from  experimental  data.  In  addition,  to  compare 
theoretical  and  experimental  results  the  solution  of  the  superposition 
that  describes  the  decay  of  the  optically  generated  carriers  is  imple¬ 
mented  in  a  program.  All  the  programs  are  written  in  Fortran  and  imple¬ 
mented  in  a  POP  11V03-L  minicomputer.  Double  precision  and  complex 
arithmetic  are  used  whe^e  needed.  A  brief  description  of  each  program  is 
included.  Listings  may  be  obtained  on  request. 

The  absorption  coefficients  in  silicon  at  a  wavelenth  of  l.O&m  are 
determined  from  an  experimental  reflection  and  transmission  coefficients 
using  equations  (1.6)  through  (1.13).  The  parameter  is  varied  until 
the  computed  value  for  R0  and  T0  matches  the  experimental  ones.  The 
coefficients  | E 2 1 ^  and  ( E 2 1 ^  can  then  be  obtained  using  equations  (1.22) 
and  (1.23).  The  aforementioned  equations  are  implemented  in  programs 
TVER. FOR  and  VER. FOR. 

The  time-dependent  modulation  of  the  CO2  laser  signal  beam  is 
theoretically  solved  in  several  steps.  The  electron-hole  density  is 
first  averaged  as  a  function  of  distance.  This  facilitates  the  calcula¬ 
tion  of  the  ambipolar  diffusivity  which  is  a  function  of  carrier  density 
and  the  hole  cross  section  at  10.59  um.  Equations  (1.41),  (1.42),  and 
(1.48)  are  used.  Subroutine  CKQ1.F0R  computes  the  coefficients  A* ,  3^ , 
and  C|<  used  in  equations  (1.41)  and  (1.42);  subroutines  AVDIFF.FOR  and 
PARM.FOR  compute  the  average  level  of  carrier  density  and  the  ambipolar 
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diffusivity  respectively.  Note  that  the  averaging  of  the  carrier  level 
is  only  done  to  obtain  average  values  for  the  ambipolar  diffusivity  and 
the  hole  cross  section.  Program  BIG. FOR  finally  calculates  the  attenua¬ 
tion  of  the  CO2  laser  probe  using  an  iterative  method  where  a  1.0  pm 
increment  is  used.  This  scheme  is  used  because  of  the  nonuniform 
spatial  distribution  of  the  carrier  density. 

In  all  the  theoretical  calculations  the  value  for  the  carrier 
lifetime  and  the  surface  recombination  velocity  are  adjusted  to  fit  the 
experimental  results.  A  lifetime  of  4.0  microseconds  and  a  surface 
recombination  velocity  of  1000  cm/second  results  in  the  best  fit. 

The  following  two  program  listings  are  Basic  programs  originally 
developed  by  Mike  Poirier  at  the  ElectroScience  Laboratory  from  The  Ohio 
State  University.  They  have  been  revised  and  updated  by  the  author. 
These  programs  are  used  in  conjunction  with  the  Tektronix  4051  micro¬ 
computer  and  7912  AD  transient  digitzer  for  data  acquisition. 
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